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Abstract 



A spin-1 model, appropriated to study the competition between bilin- 
ear (JijSiSj) and biquadratic (KijS^Sj) random interactions, both of them 
with zero mean, is investigated. The interactions are infinite-ranged and 
the replica method is employed. Within the replica-symmetric assumption, 
the system presents two phases, namely, paramagnetic and spin- glass, sep- 
arated by a continuous transition line. The stability analysis of the replica- 
symmetric solution yields, besides the usual instability associated with the 
spin-glass ordering, a new phase due to the random biquadratic couplings 
between the spins. 
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I. INTRODUCTION 



The study of disordered systems has grown very fast during the last years. Among 
these systems, spin glasses have attracted much attention. One of its main caracter- 
istic is the existence of a very rugged free-energy landscape, with many minima separated 
by high barriers. It turns out that the equilibrium state of such system becomes hardly 
accessible in an experiment, as one may guess. The spin-glass mean-field theory is well 
estabhshed 0J^, being highly nontrivial. However, the effects of fluctuations around the 
mean-field solution are very difficult to take into account in general cases, with most of 
the results been obtained for the Edwards- Anderson model [^. 

Many other spin-glass models have been investigated within the mean-field level (for 
reviews, see Refs. [^]-|^ ). Recently, much effort has been dedicated to understanding the 
phase behavior of spin-1 Ising glasses pj-^, as promising models to describe real systems 
which present multicritical phenomena. Other models which can be mapped onto spin-1 
Ising glasses were also studied recently [T0|-|T7| . However, as far as we know, none of those 
works addresses to the competition between bilinear and biquadratic random interactions. 



In fact, a few years ago a spin-glass version of the Blume-Emery-Griffiths model [IS 



was introduced in order to describe disordered magnetic lattice gases ||T9|-p3|, including 
both bihnear and biquadratic random couplings; however, due to the fact that a replica 
stability analysis was not performed, an important ingredient was missing, i.e., broken 
ergodicity, usually associated with irreversibility effects. The purpose of the present work 
is to fill this gap by investigating the overall behavior of a system which presents the 
aforementioned random interactions in a simple spin-1 model. In order to determine 
the free-energy density and the corresponding equations of state, we will use the replica 



mean-field approach. Under the replica-symmetry assumption the system exhibits 



only two phases separated by a continuous transition line. However, the stability analysis 
of the replica-symmetric solution performed within the approach proposed by de Almeida 
and Thouless |^ suggests the existence of three distinct phases. The paper is organized 
as follows. In Section II we describe the model and obtain the replica free energy. The 
replica-symmetric solution, as well as the corresponding phase diagram is investigated in 
Section HI. The stability analysis of the replica-symmetric solution is performed in Section 
IV. Our findings are summarized in Section V, where we also present our conclusions. 



II. THE MODEL AND ITS FREE-ENERGY DENSITY 

In this paper we consider an infinite-range interaction spin-glass model described by 
the Hamiltonian 

^ = - E J^As, - E K,,s^sl (1) 

where each spin Si (i = 1, 2, ■ ■ ■ , A^) can take the values —1, and 1 and the summa- 
tions are over all distinct pairs Both couplings are quenched, independent random 
variables, following probability distributions 
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where X stands for either J or K. There are two obvious hmiting cases of this modeL 
First, in the absence of biquadratic interactions {Kij = for every pair we have a 

conventional spin-1 spin-glass model. The properties of this model are quite analogous to 
those of the Sherrington-Kirkpatrick model ||5|-|9|,|2B|-P^] ; the system presents a continuous 



transition from a paramagnetic to a low-temperature spin-glass phase where the ergodicity 
is also broken [^. On the other hand, if Jij = for every pair the system becomes 



equivalent to the discrete quadrupolar-glass model investigated in Ref. [0. In this case, 
there is no sharp transition to a low-temperature phase; however, a stability analysis 
shows that in fact there is a phase transition to a low-temperature nonergodic region ||1 1|| . 
We will be most interested in the case where both J^j and Kij are distinct from zero, in 
order to appreciate the effects of their competition on the phase diagram. 
The free-energy density for this system is given by 



/5/ = -lim^, (3) 

where the bar denotes an average over the disorder. Such an average is performed by the 
so-called replica method, through the identity 

— 1 

In Z = lim , (4) 

n^O n 

which avoids the difficulty of averaging the logarithm. Using standard procedures 
we obtain 

Pf = lim -mm gn{qai3,Qai3, Pa), (5) 

where 

9n{qap, Qap,Pa) = ^ E {f^'j'dlf^ + P^^^Qlp) + ^^J' + E^' ' InTr exp(K//) 

(6) 

and 

neff = ^ E ^^^s-s^ + ^ E Q./^ isn' {s^y + ^{J' + K^) T.pI (sn' , (r) 

with the indexes a and /3 running from 1 to n. Stationarity of Qn with respect to q^^^ Qa/s 
and Pa gives the equations of state. 



P. ={is-f)^ 
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where ( denotes an average with respect to the "effective Hamiltonian" in Eq. (|^. 
Whereas the order parameters qap and Qais are already expected, the free energy depends 
also on Pa, a disorder induced order parameter which measures the fraction of spins in 
the states 5*" = ±1, for each replica a. 

In the following two sections we consider the replica-symmetric solution and its cor- 
responding stabihty analysis. 

III. REPLICA-SYMMETRIC SOLUTION 

The simplest solution of the saddle-point equations is the replica symmetric Ansatz, 
which consists in assuming 

Pa = p, Va 

qai3 = q, V(a/3) (9) 

Qa/3 = Q, V(«/5). 

Inserting this Ansatz into Eqs. (^-(0) and performing some simple Gaussian transfor- 
mations, the free-energy density becomes 



where 



and 



z(x,?/) = 1 + 2exp(A) cosh(/3jygx) , (11) 
A = ^ (p - g) + ^ (p - Q) + PK^y, (12) 



y))).y = r ^ r ^ exp(-^^)Mx, y). (13) 

" J— CO y ATX J -oo V ZTT ^ 



For the equations of state one gets. 



p =((y.2(x,t/))),,, (14) 

? =((^?(^,y))).,, (15) 
Q = {{^l{x,y)))^^, (16) 



with 



'^2(^'^) = z{x,y) cosh(/3J^x) 
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As mentioned before, the K = case is equivalent to D = in the model studied 
by Ghatak and Sherrington |2^. In this particular case it is well known that we have a 



continuous transition from a paramagnetic to a spin-glass phase at fc^T/J = 0.7901 ■ ■ ■. 
In general, the above equations present a trivial solution q = but with p 0, Q ^ 
for any temperature. In the low-temperature regime we also find another phase, with all 
order parameters distinct from zero. The phase boundary separating these two phases 
can be easily obtained from an expansion in powers of q, in either the free-energy density, 
or the equation of state for q. In either way, we find a critical frontier given by 

Q = {kBT/jf. (19) 

From the same expansions, we also ruled out the possibility of first-order transitions 
and tricritical behavior. The condition given by the above equation involves both order 
parameters Q and p, which should satisfy Eqs. ([T^) and (|16D with q set to zero. We 
are thus left with a set of three coupled nonlinear equations, which, except for some 
particular limits, has no analytical solution. We performed a detailed numerical study of 
these equations in order to check for the possibility of other types of orderings, but we 
found none, besides those already described. As a result of our analysis, we found the 
critical frontier shown in Fig. |l|. We performed an expansion for K ^ J, and verified 
that assymptotically, such a critical frontier approaches the limit ksT / J = 0.7876 ■■ ■, 
with p = Q = 0.6204 • • •. It is important to mention that within the present analysis, 
the high-temperature phase should be identified as an extension of the paramagnetic one 
already present when K = 0. This can be justified by the following argument: the free- 
energy density, as well as the order parameters p and Q, may be expanded as power series 
of K, for small values of K. Therefore, no anomalous behavior on the thermodynamical 
functions can be seen as we let K —>■ 0. Similarly, the low-temperature phase should 
be identified with the spin-glass phase ocurring at i^' = 0. Thus, as far as the replica 
symmetric solution is concerned, the biquadratic random coupling does not bring any new 
physics to this system. 

In the following section we will consider the stability of the above-mentioned solutions 
against replica fluctuations, and it will be shown that this study leads to an important 
modification on the paramagnetic side. 



IV. STABILITY OF THE REPLICA-SYMMETRIC SOLUTIONS 



Since the work of de Almeida and Thouless it is generally believed that replica- 
symmetric solutions are unstable under small fluctuations on the whole replica space. In 
our case, these fluctuations are governed by the Hessian matrix. 



G 



/ d^gn 


d^gn 




dpadpp 


dpadqu^ 


dpadQu^ 








dqu^dpa 


dqapdq^y 


dqapdQv^ 








\ dQy^dpa 


dQ^^dqafs 





(20) 



6 



where Qn is given by Eq. (j^). Stability requires that all eigenvalues of this matrix, 
evaluated within the replica-symmetric solution, should be positive (see Appendix A for 
the computation of such eigenvalues). In the limit n — we get three longitudinal 
eigenvalues, as the roots of the secular equation 



A-B- A(^) V-C T-S 

2{C-V) g - AH + ?>! - \^^^ J-A1C + ?>C 
2{J^-£) J-AIC + 3C M-4M+30- A^^) 



f211 



and two transverse ones, given by 

^^^^ g -2n+I + M-2^+0± ^[g-2n + I-M + 2Af-Of- 4[J - 2/C + C]^ 
~ 2 

(22) 

where the quantities A, . . . ,0 are defined in Appendix A. In the paramagnetic phase, 
where q = 0, one of the longitudinal eigenvalues becomes 

whereas the other two are given by 



A-B + M-4Af + 30±^[A-B-M+4Af- 30]^ - 8(^ - J^y 

•^2,3 = 2 ■ ^'^^^ 

Let us first consider the behavior of the above eigenvalues throughout the paramagnetic 
phase. Using both analytical and numerical calculations, we find that all three longitu- 
dinal eigenvalues are positive, with A^^^ vanishing along the paramagnetic to spin-glass 



transition line. However, considering the transverse eigenvalues of Eq. (|2^), we notice 
that one of them, denoted by Ai (corresponding to the plus sign before the square root), 
becomes identical to Ai^'' everywhere in the paramagnetic phase, including the critical 
frontier paramagnetic/spin-glass, where it also vanishes. Besides, in the paramagnetic 
region the second transverse eigenvalue is given by 

X^P =M-2^+0 . (25) 

Our numerical analysis shows that as the temperature decreases, and for high enough 

(T) 

values of K, A2 becomes negative, throughout the paramagnetic phase. This suggests an 
onset of irreversibility in the paramagnetic phase, associated with an ergodicity breaking, 

(T) 

as we cross the line given by A2 =0. This effect is brought about by fluctuations on the 
order parameter Qa/s, which in turn was generated by the random biquadratic couplings. 
We identify this region as a new phase, which we will call biquadratic spin-glass phase, 
with replica-symmetry breaking associated to the parameter QajS', this region should. 



then, be properly described by the Ansatz of Parisi |29|. We have also found that the 
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boundary paramagnetic/biquadratic spin-glass (where A2"^^ = 0) is a straight hne with 



slope ~ 0.077; such a numerical result is in full agreement with the one found in Ref. Wl 
{ksT/ J' f» 1.38), if one considers the proper changes of spin variables and summations in 
the Hamiltonian of Ref. [0] (which leads to K = 18 J'). 

We have also investigated numerically the behavior of all five eigenvalues in the spin- 

(T) 

glass phase. The transverse eigenvalue is negative through the whole spin-glass 
phase; this means that irreversibility is also present in this phase and so, a solution with 
replica-symmetry breaking should be employed. 

The phase diagram resulting from this analysis is shown in Fig. 0. The paramagnetic 
to spin-glass as well as the paramagnetic to biquadratic spin-glass phase boundaries should 
remain valid under a Parisi-like treatment. However, it is possible that the biquadratic 
spin-glass to spin-glass frontier changes under replica-symmetry breaking in both matrices 
Q and q. Thus, the corresponding boundary shown in Fig. 0, which was obtained within 
the replica-symmetric solution, should be seen as a rather schematic one, although there 
is no physical reason to expect a substantial qualitative change. The correct treatment 
based on Parisi's Ansatz is very difficult in this case, since it involves nonlinear integro- 
differential equations at finite temperatures, which are hard to solve numerically. Such 
an analysis is beyond the purpose of this paper. 



V. CONCLUSIONS 

We have studied a solvable spin-1 model, including both bilinear and biquadratic 
random exchanges, with zero means and variances J and K, respectively. The model was 
solved through the replica formalism. Three types of order parameters were introduced 
to describe the system in the replica space: a density (pa), which measures the fraction of 
spins in the states S"" = ±1 and two spin-glass-like matrices, represented by the bilinear 
and biquadratic matrix elements qa/s and Qa/s, respectively. 

The replica-symmetric solution leads to a continuous transition from a high- 
temperature paramagnetic phase to a low-temperature spin-glass phase, signaled by the 
onset of the spin-glass order parameter q. The corresponding critical frontier is almost 
temperature-independent, especially for large values of the variance K. We have also 
analysed the eigenvalues of the stabihty matrix associated with fluctuations around the 
replica-symmetric solutions. We verified numerically that one of the replicon eigenvalues 
is always negative throughout the whole spin-glass phase, implying an instability of the 
replica-symmetric solution. We have also noticed that the paramagnetic phase presents a 
similar instability (associated with the matrix elements Qa/s) , for sufficiently large values 
of the variance K, giving rise to a new phase which we have called biquadratic spin- 
glass phase. Such instabilities may be related to the onset of irreversibility effects, i.e., 
the response functions could depend on the history of the system (e.g., field-cooling and 
zero-field-cooling measurements may lead to different results) 0,^. On the basis of our 
findings, many of them from numerical analysis, we conclude that for K ^ the system 
presents at least three distinct phases, in which two of them should be properly described 
through a replica-symmetry breaking procedure. The frontier separating the biquadratic 
spin-glass and spin-glass phases requires further investigation. It is not clear if a full Parisi 
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solution would change substantially its location. We hope to address to this point in a 
future work. 
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APPENDIX A: STABILITY ANALYSIS OF THE REPLICA-SYMMETRIC 

SOLUTION 



The elements of the Hessian matrix G, defined in Eq. (^), are given by 
dpadpf} 



d'gn 

dqapdq^s 

dqafsdQ-ys 

d^gn 
dQafsdQrys 

where 

For the rephca-symmetric solution the eigenvectors of G have the form 



u 



where 



— Pa — P, Val3 — QajS — <?, Ca/B — QajJ — Q, 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 

(AT) 



(A^ 



(A9) 



represent Gaussian fiuctuations. Following de Almeida and Thouless pH|, we start with 
the eigenvector totally symmetric under replica-index permutations 



ea = a, r]ap = b, = c, for a, /? = 1 . . . n, 



(AlO) 



which correspond to the longitudinal eigenvectors, according to the conventional classifi- 
cation |3^. For a finite value of n, the corresponding eigenvalues follow from 



A^^^a = Aa+{n- l)Ba + {n - l)Cb + ^(n - 2)(n - l)Vb 

-f(n - l)8c + -{n -2){n- l)Tc, 
A(^)6 = 2Ca + {n- 2)Va + Qb + 2{n - 2)m + (^"^K^-^) -^^ 



(All) 



(A12) 
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where 



+ + 2(n - 2)/Cc + 3) ^^^ 
A(-^)c = 2£:a + (n - 2):ra + J^fe + 2(n - 2)/C6 + (^-2)(^-3) ^^ 
+Mc + 2{n - 2)Mc + (^-2)(n-3) ^^^ 



C 

£: 
c 

M 
N 
O 



d^gn 



dpadpa 



dpadpp 
d'gn 



RS 



RS 



1 - 



2 



{i-p)p 



- Q) , 



dpadqafi 



dpadqp^ 



RS 



RS 



{(5Jf{(5Kf 



(pq-w), 



dpadQa(3 



dpadQp^ 



dqapdqap 



dqapdqaj 
d'gn 



dqapdq-ys 
d'gn 



RS 



RS 



RS 



RS 



RS 



(p-i)Q, 



{pQ-W), 



dqapdQap 



dqapdQ 
d^gn 



dqapdQjs 
d'gn 



RS 



RS 



RS 



dQapdQafi 

d'^gn 



dQapdQ 

oc~i 

d^gn 



dQafidQjS 



RS 



RS 



RS 



{(5Jf [l + {(5J)\q'-Q) 
iPJnq'-w), 

-{PJ)\PK)\Qq-v), 

= {PKf [l + {PKf{Q-l)Q 

= {(3Kf (g^ - w) , 
= (PK)' {q' - S) . 



(A13) 



(A14) 
(A15) 
(A16) 
(A17) 
(A18) 
(A19) 
(A20) 
(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 
(A27) 
(A28) 



In the above equations, different rephca labels are distinct from one another, whereas ..\rs 
means elements of the Hessian matrix evaluated within the replica-symmetric solution; 
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the parameters p, q and Q are given by Eqs. (|T^-(p!6D and 



S={{^t))^^, (A29) 

v = {{^l^l))^^, (A31) 

W = {{^1))^^, (A32) 

w = {{^l^2))^. (A33) 



The next step is to find the anomalous eigenvalues. These correspond to breaking the 
symmetry of the vector u with respect to one specific replica-index, denoted herein by 6, 




1]°'^ = hi for a or (3 = 6 

r]"(^ = b2 for a,P^e 

= Ci for a or P = 9 

= C2 for a,p^e 



(A35) 



(A36) 



Orthogonality with the replica-symmetric eigenvector implies 

a2 = r, &2 = C2 = -. (A37) 

n — 1 n — 2 n — 2 

Thus, the anomalous eigenvalues follow from 

A(^)ai = {A- B)ai + (n - 1)(C - V)bi + [n - 1){S - ^)ci, (A38) 

A(^)6i = ^^^—^{C - V)ai + [g + nn-{n- 3)J]6i 
n — I 

+ [J + nlC- {n-3)C]ci = 0, (A39) 
A(^)ci = - + [J + n}C-in- 3)C]bi (A40) 

+ [A^ + riAT - (n - 3)0]ci = 0. 

In the limit n the longitudinal and anomalous eigenvalues coincide and may be 
obtained from Eq. (^I]) . 

Finally, we must find the transverse eigenvalues. In this case, two replica indices are 
fixed and the symmetry between replicas is broken with respect to such indices (herein 
denoted by 9 and z/). One has, 

{e" = as for a = 6^ or 1/ { \a-\\ 

= 04 for a^e,iy ^ ' 
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' rj°'f^ = 63 for a or /9 = ^ or 1/ 

< T]^'' = T]""" = 64 for a^e,iy (A42) 
r]'^P^h^ for a,P^e,iy 



' = C3 for a or P = 9 or u 

< ^^" = ^^° = C4 for a^e,p (A43) 
^"/? = C5 for a,P^e,u 

These eigenvectors must be orthogonal to both longitudinal and anomalous ones. Thus, 
it follows that 03 — 0,4. = and 

63 = -{n - 2)64= ^{n - 2){n - 3)h, 

C3 = -(n-2)c4= ^(n-2)(n-3)c5. (A44) 
Prom these observations, the transverse eigenvalues are given by the secular equations 

X^^^h = {Q-2n + I)h + {J-2IC + £)C3, (A45) 
A(^)c3 = ( J - 2/C + C)h + {M-2Af + 0)c3, (A46) 

which are independent of n. 
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FIGURES 



FIG. 1. Phase diagram obtained within the replica- symmetric solution, presenting a contin- 
uous transition from a high-temperature paramagnetic phase (P) to a low-temperature spin-glass 
phase (SG). 



FIG. 2. Phase diagram resulting from the stability analysis of the replica- symmetric solu- 
tion. This solution is stable throughout the paramagnetic (P) phase only. The biquadratic 

spin-glass phase (BSG), which occurs for large values of K , is characterized by an instabil- 
ity of the replica-symmetric paramagnetic solution. In the whole spin-glass phase (SG) the 
replica- symmetric solution is unstable. 
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